"""使用保存下来的数据进行一些数据处理和绘制工作"""
import matplotlib.pyplot as plt
import scipy.io as scio

# 绘制21阶，正弦\Delta C的仿真结果，以确定仿真时的相对容差和绝对容差
# Ip1和Ip2~1e-4, Is~1e-8, Qeq1~1e-11, Qeq2~4e-18
# Uoa和Uob~1e-3, dUoa/dt和duob/dt~1.5e4, Ulp1e-2, dUlp/dt~5e-3
# 参数是按Mance论文来的，补偿电容阵列5个但都未启用
M_results = scio.loadmat(r'D:\WorkSpace\temp\data\trying2.mat')
M_t = M_results['T']
M_t.shape = M_t.size,
M_Y = M_results['Y']
# Is比其他两个都小很多，为了看清分开画
fig, axes = plt.subplots(nrows=2, ncols=3)
axes[0,0].plot(M_t, M_Y[:, 0],label='Ip1')
axes[0,0].legend(fontsize=8)
axes[0,1].plot(M_t, M_Y[:, 1],label='Ip2')
axes[0,1].legend(fontsize=8)
axes[0,2].plot(M_t, M_Y[:, 2],label='Is')
axes[0,2].legend(fontsize=8)
axes[1,0].plot(M_t, M_Y[:, 3],label='Qeq1')# 这个仿真里没有开补偿电容
axes[1,0].legend(fontsize=8)
axes[1,1].plot(M_t, M_Y[:, 4],label='Qeq2')# 可以认为Qpri1和Qpri2的量级和这个类似
axes[1,1].legend(fontsize=8)
# 
fig, axes = plt.subplots(nrows=2, ncols=2)
axes[0,0].plot(M_t, M_Y[:, -2],label='Ulp')
axes[0,0].legend(fontsize=8)
axes[0,1].plot(M_t, M_Y[:, -1],label='dUlp/dt')
axes[0,1].legend(fontsize=8)
axes[1,0].plot(M_t, M_Y[:, -6],label='Uoa')
axes[1,0].plot(M_t, M_Y[:, -4],label='Uob')
axes[1,0].legend(fontsize=8)
axes[1,1].plot(M_t, M_Y[:, -5],label='dUoa/dt')
axes[1,1].plot(M_t, M_Y[:, -3],label='dUob/dt')
axes[1,1].legend(fontsize=8)
plt.show()
